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ABSTRACT 

This  study  presents  the  derivation  of  a  diffuser  design  parameter, 
which  takes  into  account  losses  due  to  the  geometry  of  the  diffuser. 

A  criterion  for  design  of  annular  diffusers,  using  the  design 
parameter  to  predict  losses,  is  developed  from  reference  data  and 
model  tests  . 

A  curved  wall  annular  diffuser  was  designed  and  tested  to  verify 
the  predictions  of  the  design  criterion  for  an  arbitrary  diffuser. 
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I.   INTRODUCTION 

In  the  design  of  annular  diff users,  such  as  those  found  in  jet 
engines,  and  turbomachines ,  there  has  been  a  problem  defining  param- 
eters which  indicate  sufficiently  the  optimum  design  range  for  varying 
geometries.   Many  schemes  have  been  devised,  such  as  the  so-called 
"equivalent  conical  angle  method",  to  relate  a  diff user  to  a  known 
shape.   The  difficulty  in  these  approaches  has  been  that  each  person 
using  the  method  developed  a  slightly  different  approach  which  resulted 
in  a  variety  of  descriptions  for  the  same  diffuser.   Gleason  [1]  gives 
a  sunKBS.ry  of  such  methods  and  points  out  the  difficulty  in  correlating 
data  from  different  sources. 

Some  apparent  problems  arise  when  describing  a  diffuser  by  inlet 
area,  exit  area,  length,  difference  in  inlet  radii,  equivalent  cone 
angle,  or  combinations  of  these  parameters  developed  by  empirical  means 
Although  these  terms  give  some  intuitive  feel  for  the  way  a  diffuser 
would  perform,  two  diffusers  which  differ  widely  in  performance  could 
have  the   same  definition. 

Tfae  simple  example  of  Figure  1,  where  the  two  diffusers  have  the 
same  Inlet  area,  exit  area,  length  and  possibly  the  same  divergence 
angle,  makes  it  immediately  obvious  that  the  two  diffusers  would  not 
have  the  same  losses.   The  velocity  and  pressure  distributions  through 
the  diffuser  channel  have  a  direct  effect  on  the  boundary  layer 
development,  which  in  turn  controls  the  frictional  losses.   The 
widely  differing  geometries  of  the  two  above-mentioned  diffusers 
would  cause  large  differences  in  boundary  layer  growth.   Therefore, 
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FIGURE    1 

Example    of  diffusers  which   previously  had   the 
same   description. 
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the  performance  of  the  two  would  differ  and  a  similar  description 
would  not  be  appropriate. 

Sovran  and  Klomp  [2]  have  empirically  developed  a  fairly  successful 
set  of  parameters  for  straight  wall  annular  diffusers.   The  geometric 
area  ratio  was  taken  as  the  primary  factor  affecting  diffusion,  and  a 
non-dimensionalized  length  when  combined  with  the  area  ratio,  gave  an 
overall  criterion  for  the  pressure-gradient  affecting  the  boundary 
layer  development.   Their  description  of  annular  diffusers  is  at 
present  limited  to  designs  with  straight  walls. 

Vavra  [3  and  A]  presented  the  idea  of  a  shape  factor  or  parameter 
which  includes  the  effects  of  diffusion  and  boundary  layer  influence 
for  the  evaluation  and  preliminary  design  of  diffusers  of  arbitrary 
shapes . 

The  objectives  of  this  study  were  to  present  the  derivation  of  the 
shape  parameter,  Q,  to  substantiate  its  usefulness  from  reference  data, 
and  to  evaluate  the  validity  of  the  results  by  construction  and  investi- 
gation of  a  family  of  straight  wall  diffusers  such  as  those  in  Ref.  2. 
To  show  that  the  design  shape  parameter  may  be  applied  to  diffusers 
other  than  those  with  straight  walls,  a  curved  wall  diffuser  such  as 
might  be  found  at  the  exit  of  the  high  pressure  compressor  of  a  jet 
engine  was  designed  for  an  optimum  Q  and  evaluated  by  tests. 
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II.   DEVELOPMENT  OF  THE  SHAPE  PARAMETER  OMEGA  (Q) 

Tfae  function  of  a  diffuser  is  to  convert  a  portion  of  the  kinetic 
energy  of  the  free  stream  into  static  pressure  by  a  compression  process 
similar  to  that  in  a  mechanical  compressor.   The  shape  parameter, 
Omega  (Q) ,  to  be  used  for  diffuser  design,  can  be  evolved  by  examining 
the  differential  entropy  changes  of  such  a  polytropic  compression 
process  . 

A.   ENTROPY  CHANGE  OF  A  ONE -DIMENSIONAL,  FRICTIONAL  FLOW 

An  -.expression  for  the  differential  entropy  change  of  a  polytropic 

compression  may  be  found  by  investigating  a  one -dimensional  flow  in  a 

diverging  channel,  such  as  shown  in  Fig.  2.   A  momentum  balance  on  the 

differential  volume  of  Fig.  2  gives 

m  (V  /  dV)  -  mV  =  pA  -  (p  /  dp)  (A  /  dA)  /   (p  /  dp/2)  dA  -  TCdL 

/  m 

which,  reduces  to  /^\  ^  ' 

m  dV  =  -  Adp  -  TCdL  (2) 
or 

dp  =  -  f  dV  -  &  dL  (3) 
where    the  values    of  A,   C,    p,   dL,   V,    and  T   are   as  defined   in  Fig.    2. 
Introducing   the   equation   of   continuity, 

m  =  pVA  (4) 
into   equation    (3)    results    in 

dp  =   -   pVdV   -  ^  (5) 


or 


f  =  -  d(V2/2)   -  I  £  dL  (6) 
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FIGURE  2 
One -dimensional  flow  in  a  diverging  channel 
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with  the  First  Law  of  Thermodynamics 


or 


dq  =  du  /  pdv  =  dh  -  vdp  =  dh  -  -  dp  =  Tds    (7) 


^  =  dh  -  Tds  (8) 

P 

Equating  (6)  and  (7)  gives 

dh  /  d(V2/2)  -  Tds  =  -  7  dL  (9) 

p  A 

The  first  two  terms  of  (9)  combine  to  the  derivative  of  the  total 
enthalpy,  H,  and  if  the  flow  is  restricted  to  an  adiabatic  case,  with 
no  energy  in  any  form  transferred  between  the  fluid  and  its  surroundings, 
the  total  enthalpy  is  constant  and  its  derivative  is  equal  to  zero. 

Therefore, 

dH  =  d  (h  /  V2/2)  =  0 
and  equation  (8)  simplifies  to 

Tds  =  1  7  dL  (10) 

P  A 

Introducing  the  concept  of  the  local  coefficient  of  skin  friction 

in  a  one -dimensional  flow, 

cf  = T-J  (11) 

f    l/2Pr 

equation  (10)  becomes 

c 

Tds  -  y-   V2  |  dL  (12) 

Introducing  equation  (A)  and  the  equation  of  state  for  a  perfect  gas 

in  (12)  results  in  an  expression  for  the  differential  entropy  change, 

c.    .2  c_    .2   R 

,  fm  C,T  fmgC  /ion 

ds  =r~v  T-dL  =  r^r;3 dL  (13) 

p  A  r     A 
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B.   ENTROPY  CHANGE  FROM  THERMODYNAMIC  CONSIDERATIONS 

The  differential  change  in  entropy  can  also  be  found  for  an 
adiabatic,  polytropic  compression  process  as  shown  in  the  T-s  diagram 
of  Fig.  3.   The  definition  of  polytropic  efficiency  is, 

dT. 

The  isentropic  relation  for  a  differential  compression  from  p  to 

p  +  dp  is 

2-1 

T  +  dT.    r     .    ,  N  y 

— t"^  ■  C^T2)  <13> 

By  expanding  (dp/p)  by  a  binomial  series,  and  assuming  dp  to  be  small 

compared  to  p, 

dT.   =T^^  (16) 

is      7   p 

Combining  (14)  and  (16)  gives 

dT  „  I_  2Z1  §£  (17) 


1p  7 


Solving  equation  (7)  for 


and  introducing 


A  dT     l        A 

ds  =  c  —  -  — —  dp 
p  T    pT   p 


c  -  R  -hr 
P    g  7-1 


and  tihe  perfect  gas  law 

p/p  =  Rg  T 
results  in  an  expression  for 


g  1-7- 1  T 
Substituting  (17)  into  (18)  gives 


--I      fef-?]  (18) 


ds  =  r  &  [jl  .  {\  (19) 

g  P    Ltj, 


p 
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FIGURE  3 

T-s  diagram  for  an  adiabatic,  polytropic  compression  in 
a  stationary  duct. 
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C.   DEFINITION  OF  OMEGA  (Q) 

Equating  the  two  values  of  differential  entropy  change  from 
equations  (13)  and  (19)  results  in 


c.  .  _  R       R  ri 

2    p  A3  P  P    Lr,p 


dp 


(20) 


Introducing  the  concept  of  referred  flow  rate,  a  dimensionless  quantity 

=  v  g   1   • 


Ir    PlAj_ 


and  substituting  (21)  into  (20)  and  rearranging  gives 


(21) 


For  a  polytropic  process 

p/p  =  constant 
where  n  is  the  polytropic  exponent. 
Rewriting  (22)  in  the  form 


1_ 


pi     Vpr 


■) 


(22) 


2"mlr 


<£>  h  «.  -  fr  - 1]  &T-  £)  C") 


1      -"p    ^  T 
gives  two  independent  differentials  which  can  be  integrated 

Defining  the  left  side  of  (23)  as  dX..  and  integrating, 


2  A\ 


i  - jr^ivQr)  k[dL 


(24) 


where  L  is  the  meridional  length  along  the  mean  streamline. 

Assuming  an  average  value  of  coefficient  of  friction,  cf,  for  the 
whole  surface  and  that  m.   is  a  constant  for  a  given  flow  rate,  (24)  may 
be  rewritten  as 
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L     A 
Xl   =  2~  ™2    ^^ 


lr  J  \A  J     A 
0 


(25) 


The    integral   of    (25)    is   a   parameter  which   contains    the    flow  area   ratio, 
the   diffuser    length,    and   the  wetted   surface,    and  will   be    called   the 
"diffuser   shape    parameter",    defined  by 


0 


o 


dL 


(26) 


D.      RELATIONS   BETWEEN  Q  AND   DIFFUSER  PERFORMANCE 


Defining   the    right-side   of    (23)    as   dX    , 


p2/pl 


X, 


P2/Pl  ri  -,    .       .1/n         .    x 


If  n      and  n  are   assumed   constant, 
'P 


X2   "  [J"  "    ll    W^  [(^)(1/"  +   1J-    0        <28> 


For 


p2   =  Pl  +  Ap 
the  binominal  expansion 


-  1    + 


1/n  +   1 


1  + 


G*0 


^  . 


Assuming  Ap  is  much  less  than  p.  and  omitting  higher  order  terms  from 
the  binomial  expansion,  equation  (28)  becomes 


-E 


^£ 


2   <-ri     J  P 


(29) 
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Equating  X,  and  X  , 


& 


Ap    f  •  2  n 
j  ?l        2        lr 


With  the  equation  of  continuity,  equation  (21)  becomes 


K 


lr    111  p..A 
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Substituting  m,   into  equation  (30) 


a 


_Ap_ 


1/2d1V1 


2  =  c  Q 


For  a  constant  total  temperature,  equation  (18)  becomes 

ds  =  -  R 


8  P* 


(30) 


(31) 


(32) 


Thus  equation  (19)  becomes 

dP. 


Integrating  (33) 


In 


QO-fr^]-© 


(33) 


(34) 


Using  the  assumption  that  P   =  P   +  AP  and  p_  =  p  +  Ap,  where  AP 

t2    tl 
and  Ap  are  small  compared  to  P  and  p,  respectively,  and  with  the 

truncated  power  series  expansion  of  In  (1  +  x)  =  x,  (34)  becomes 


AP. 


& 


-  1 


Ap 


J  p. 


or 


[i-  -  l]  Ap  =  APt  (£-)  =  4P 

'P  tl 


1  +  2  V 


(35) 
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2 

Since  for  an  incompressible  flow,  M  ,  is  small  compared  to  unity, 


[i-  -  l]  Ap  -  AP£  (36) 


P 
Substituting  (36)  into  (32)  gives 

APt 

Cfn= — j  (37) 

f     1/2  Pl  VLZ 

In   Ref.    2,    the   Pressure  Recovery  Coefficient    is   defined  as   the   ratio 
of   the   actual   pressure   rise   and   the   maximum  attainable    one    for   the 
particular    flow  rate,    or 

__     P2-PX 


'pr        -.-       ..2 


(38) 


1/2   PjVj 

Defining  the  Ideal  Pressure  Recovery  Coefficient  as  that  which 
would  result  for  an  isentropic  compression  producing  the  theoretical 

pressure  p9   of  Fig.  3, 
th 

P2th"Pl 
C   .  «  — ^ 5"  (39) 

pri    1/2  plVl2 

For  isentropic  flow  Ap  =0,  and  assuming  density  to  be  constant  and 
introducing  the  equation  of  continuity, 

p/2  V  2  -  p/2  V  *       Ji^J 

c   . 1 2-  =  1  -  (-i)         (40) 

Pri       p/2  Vx2  ^2-" 

From  equation  (40), 

p/2   V22   =  p/2   VL2    [1   -  Cpr.] 

Then 

P  _P   =  p  +  p/2  VL2  -  p2  -  p/2  VL2  [1-C    ] 
Cl  C2  V 

or 
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p.  -  p. 

?        Pt    "   P? 

+  C      .  (41) 


ll        H        pl    "    p2 


p/2   V*        p/2   Vx2  pri 

Using   the   definition  of  C      ,    (38),    and  equation    (37)    there    is 

dc\  =  C      .    -   C  (42) 

f  pri  pr  v      ' 

or 

C      .    -   C 

\  -  ^r^  («> 

Hence,  after  having  previously  determined  Q,  the  average  frictional 
loss  coefficient,  c  ,  can  be  calculated  from  the  measured  quantities 

pl'  ql'  Al'  P2'  and  A2' 
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III.   TEST  INSTALLATION  AND  INSTRUMENTATION 

The  apparatus  used  in  this  investigation  was  located  in  Building 
230  at  the  Turbopropulsion  Laboratories  of  the  Department  of  Aeronautics, 
Naval  Postgraduate  School,  Monterey,  California. 

A.   AIR  DELIVERY  SYSTEM 

A  schematic  drawing  of  the  installation  is  shown  in  Fig.  4.   The 
air  source  was  a  200  H.P.',  Carrier  centrifugal  compressor  with  a  maxi- 
mum flow  rate  of  about  four  pounds  mass  of  air  per  second  at  a  pressure 
ratio  of  about  two.   Gross  throttling  was  accomplished  by  a  main  by-pass 
valve  ahead  of  the  air  coolers.   Fine  adjustments  of  the  flow  rate  were 
made  by  bleeding  air  from  plenum  (1)  through  a  secondary  by-pass  valve. 

The  air  from  plenum  (1)  discharged  through  a  flow  straightener  into 
a  four  inch  pipe  containing  a  flow  measuring  nozzle.   The  nozzle  used 
was  the  same  as  that  described  by  Kelly  [5],  which  had  a  diameter  of 
3.022  inches  and  a  diameter  ratio,  3,  of  0.734. 

The  flow  was  diffused  into  plenum  (2)  by  a  ten  degree  conical 
diffuser  extending  downward  from  the  inlet.   The  discharge  of  plenum 
(2)  was  into  a  rectangular  wooden  duct  containing  flow  straighteners 
and  a  total  temperature,  iron-constantan  thermocouple  probe. 

The  mean  stream  velocity  in  the  duct  at  maximum  flow  rate  was 
approximately  20  to  25  feet  per  second.   Transition  pieces  shown  in 
Fig.  5  were  used  to  change  the  geometry  of  the  flow  channel  from 
rectangular  to  circular  and  to  accelerate  the  flow  to  the  diffuser 
inlet  velocity  and  pressure. 
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B.   DESCRIPTION  OF  DIFFUSER  MODELS  AND  ASSEMBLY 

The  diffuser  assemblies  of  the  two  models  constructed  and  their 
supporting  assembly  are  shown  in  the  photographs  of  Figures  6  through 
10.   The  supporting  assembly  consisted  of  the  duct  transition  pieces 
with  a  cylindrical  portion,  shown  in  Fig.  6,  containing  three  equally 
spaced  support  struts  for  a  cantilevered  center-body,  so  that  no 
obstructions  or  supports  would  be  necessary  in  the  diffuser  flow 
channel.   The  flow  was  accelerated  to  the  diffuser  entrance  in  an 
annular  space  created  by  the  nose  of  the  center-body  assembly  and  a 
converging  outer  wall. 

A  straight  section,  where  the  annular  area  remained  constant,  was 
provided  just  downstream  of  the  inlet  to  provide  room  for  flow  measure- 
ments and  surveys,  as  shown  in  the  schematic  of  Fig.  11.   Smooth 
transitions  to  the  diffuser  contour  were  arranged  as  shown  in  Figures 
8  and  9.   All  computations  of  the  diffuser  shape  parameter  included  the 
flow  area  from  the  inlet  plane  to  the  exit  plane  of  the  diffuser, 
inclusive  of  the  cylindrical  inlet  section. 

The  diffuser  inlet  plane  was  provided  with  three  supports  for 
probes,  equally  spaced  around  the  outer  periphery  and  bisecting  the 
angles  formed  by  the  supporting  struts  of  the  center-body.   These 
supports  could  accommodate  both  hot-wire  and  pressure  probes  for 
radial  surveys.   One  support  was  used  for  a  Kiel  probe  for  the  measuring 
of  total  pressure  at  the  inlet  plane.   Three  static  pressure  taps  were 
provided  in  the  same  plane,  120  degrees  apart  on  the  outer  wall.   Three 
sets  of  static  taps  were  placed  axially  downstream  from  the  inlet  static 
taps,  at  two- inch  intervals  on  the  outer  wall  of  the  first  model,  and 
on  both  the  inner  and  outer  walls  at  one-inch  intervals  on  the  second 
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FIGURE  5 
Transition  ducting,  looking  downstream, 


FIGURE  6 
Support  assembly,  including  support  struts 
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FIGURE  7 

Diffuser  support  assembly  and 
converging  outer  wall  (with 
Model  (1)  attached). 


FIGURE  8 

Diffuser  support  assembly  with 
Model  (1)  cantilevered  center-body 
installed . 
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FIGURE  9 
Center-body  of  Model  (2) 
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FIGURE  10 
Outer  wall  of  Model  (2)  with  Support  Assembly, 
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model.   Static  pressure  readings  were  found  to  be  equal  at  all  three 
stations  around  the  circumference,  and  the  three  taps  at  each  axial 
location  were  therefore  connected  together. 

Model  (1)  had  straight  walls  as  illustrated  in  Fig.  12,  with  a  wall 
divergence  angle  of  15  degrees  for  both  the  inner  and  outer  contours. 
The  construction  was  such  that  a  series  of  four  diffusers  with  differ- 
ent shape  parameters  could  be  produced  by  simply  cutting  off  the  end  of 
the  model.   Model  (2)  had  a  more  complicated  shape,  and  was  designed 
as  described  in  Appendix  B  to  show  the  validity  of  the  Q  parameter  for 
curved-wall  diffusers.   A  summary  of  the  geometrical  characteristics 
of  both  models  is  given  in  Table  I. 

C.   INSTRUMENTATION 

Pressure  readings  for  the  flow  nozzle  and  the  diffuser  inlet  were 
made  with  a  Texas  Instrument  Fused  Quartz  Pressure  Gage,  which  is 
shown  in  the  overall  view  of  the  experimental  apparatus  of  Fig.  13. 
A  constant  pressure  source  calibration  of  the  Pressure  Gage,  to  a 
Meriam  Vernier  Manometer  filled  with  water,  gave  a  constant  value  of 
0.800  inches  of  water  per  count,  with  a  reading  accuracy  of  +  0.001 
counts.   This  result  was  identical  to  that  of  Beck  [6]. 

Pressure  readings  for  five -hole  probe  measurements  were  made  on  a 
235  cm,  eight  tube  water  manometer  bank,  which  can  be  read  to  an 
accuracy  of  +  0.1  cm.   Static  pressure  measurements  along  the  diffuser 
length  were  made  on  a  92  inch,  20  tube  water  manometer  board  with  a 
reading  accuracy  of  +  0.05  inch.   All  pressure  taps  were  connected  to 
their  respective  reading  devices  with  plastic  tubing.   A  standard 
mercury  wall-barometer  was  used  for  atmospheric  pressure  readings. 
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Temperature  readings  from  the  thermocouples  upstream  of  the  flow 
nozzle,  and  in  the  wooden  duct,  were  obtained  with  a  Leeds  and  Northrup 
portable  potentiometer  with  compensated  cold- junction,  giving  a  read- 
ing accuracy  of  +  0.01  mV.   Adiabatic  flow  was  assumed  from  the  area 
of  temperature  measurement  to  that  of  pressure  measurement. 

Hot-wire  measurements  were  conducted  with  a  two  channel,  constant 
temperature,  hot-wire  anomometer  manufactured  by  Applied  Science, 
Carmel,  California.   The  wires  were  0.00015  inch  diameter  tungsten  with 
a  0.1  inch  gap.   Surveys  were  taken  radially  from  the  outer  to  the 
inner  wall. 
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IV.   EXPERIMENTAL  PROCEDURES  AND  COMPUTATIONS 

All  reduction  of  data  was  carried  out  by  the  Fortran  IV  computer 
program  "Omega"  described  in  Appendix  C. 

A.   PRESSURE 

Pressures  obtained  from  the  Texas  Instrument  Fused  Quartz  Pressure 
Gage  were  reduced  to  pounds  per  square  inch  by 

p  =  0.03613  K  pTI  (44) 

where  K  =  0.8,  the  conversion  factor  from  counts  to  inches  of  water, 
p   is  the  pressure  read  in  counts  on  the  Pressure  Gage,  and  the 
constant  value,  0.03613,  is  the  conversion  factor  for  inches  of  water 
to  pounds  per  square  inch  from  Ref.  7. 

Barometric  pressure  was  converted  from  inches  of  mercury  to  pounds 
per  square  inch  by  equation  (3-2)  of  Ref.  6. 

0.000163''!  t 


PA  =  0.4912  p 


1  "  1  +  0.0001818  t  J        (   ' 


c 


where  the  coefficients  of  t  ,  the  room  temperature  in  degrees  centi- 
grade, are  the  thermal  expansion  coefficients  of  mercury  and  brass. 
The  factor  0.4912  is  the  conversion  factor  for  inches  of  mercury  at 
0  C  to  pounds  per  square  inch  from  Ref.  7. 

B.   TEMPERATURE 

Thermocouple  readings,  recorded  in  millivolts,  were  converted  to 
degrees  Fahrenheit  by  the  temperature  relaL ionship, 

t  =  32.144  /  35.77  mV  -  0.4518  mV  (46) 
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obtained  from  a  least  squares  fit  of  data  from  Ref.  8.   This  relation 
holds  for  temperatures  below  100  degrees  Fahrenheit,  with  the  cold- 
junction  maintained  at  32  degrees,  or  compensated  to  32  degrees  by  the 
standard  cell  of  the  Leeds  and  Northrup  measuring  apparatus. 

C.   FLOW  RATE 

After  setting  the  desired  flow  rate,  the  fluid  temperatures  were 

allowed  to  stabilize  before  data  were  taken.   Upstream  static  pressure, 

p  ,  differential  pressure  across  the  flow  nozzle,  (p  -p.)>  and  the 
ru  r  u  d 

upstream  total  temperature,  tf    ,  are  necessary  to  compute  flow  rate 
as  described  in  Appendix  A  of  Ref.  5. 


p 
*   u    c 


w  =  w  /Tjr 

V  fu  V  I 
and 


(47) 


*  p„  "  ?<\ 

w  -  11.30  e  f   — (48) 

t  v   p 

•     u 

*  2 

where  w  is  the  equivalent  flow  rate ,  which  has  dimension  in  . 

The  so-called  expansion  coefficient,  c,    is 


,-  2.433  i^ 


1.428  ^,    ,.   N0.286 
1  -  (1-x) 


i  -  .aosa-x)1-42^ 


(49) 


where 


x 


pu  (50) 


and 


ft  ■  1  '  °-0015  C^IoT-)  (51) 


is  a  correction  factor  that  takes  account  of  the  thermal  expansion  of 
the  nozzle  throat  diameter. 
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D.   CHECK  OF  NOZZLE  FLOW  RATE  WITH  SURVEY  DATA 

Preliminary  surveys  at  the  inlet  plane  with  a  hot-wire  probe 
indicated  that  the  velocity  profiles  at  the  three  stations  surveyed 
were  similar.   Figure  14  is  a  least-squares  fit  of  the  data  points  from 
all  three  survey  stations.   The  data  were  taken  at  varying  flow  rates 
at  each  station,  and  the  velocities  were  found  to  vary  within  +  1 
per  cent  of  the  mean  obtained  with  the  least-square  fit.   Hence,  the 
flow  was  very  nearly  axisymmetric ,  and  the  non-dimensional  profiles 
remained  unchanged  at  the  different  flow  rates  investigated.   The 
boundary  layer  was  very  thin  at  the  inlet.   The  actual  thicknesses  were 
0.025  inch  on  the  inner  wall  and  0.75  inch  on  the  outer  wall  according 
to  the  measurements. 

To  verify  the  pressure  and  temperature  measurements,  a  comparison 
was  made  of  the  flow  obtained  from  the  survey  data  at  the  inlet  and 
that  measured  by  the  flow  nozzle.  A  survey  with  a  five-hole  probe  was 
made  at  the  inlet  plane  to  obtain  the  actual  values  of  the  velocities 
which  are  plotted  in  Fig.  15.   Data  for  the  five-hole  probe  were  cor- 
rected by  use  of  the  calibration  curves  supplied  by  the  manufacturer, 
United  Sensor  and  Control  Corporation.   The  velocity  profile  obtained 
was  essentially  the  same  as  that  of  the  hot-wire  surveys.   The  hot-wire 
values  were  used  for  the  boundary  layer  regions,  where  the  immersion 
corrections  for  the  pressure  probe  were  unreliable.   The  yaw  angles  of 
the  flow  were  less  than  +  1/2  ,  and  the  flow  had  pitch  angles  between 
three  and  six  degrees,  with  the  velocities  pointing  outward  from  the 
center.   The  greatest  velocities  occurred  at  the  outer  wall. 
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The  five-hole  probe  survey  showed  the  total  pressures  and  the 
velocities  to  be  within  two  per  cent  of  the  mean  values  in  the  turbu- 
lent central  region  outside  the  boundary  layers.   An  average  Mach 
number  of  the  flow  in  this  central  region  was  determined  from  the 

isentropic  relationship 

7-1  1/2 


The  average  static  temperature  at  the  inlet  was  determined  by 
assuming  adiabatic  flow  from  the  upstream  total  temperature  probe  to 
the  inlet  plane,  which  gave  the  isentropic  relationship 

1 
Hence , 


T 

j'1  =  1  +  Zji  «\  (53) 


Tl  =  FT"!  <54> 

1  +  •Z— =  M 
2    1 

Static  pressure  was  found  to  be  nearly  uniform  across  the  inlet 

plane,  so  that  density  was  almost  constant  at  this  station.   From  the 

perfect  gas  relationship, 

P]L(144) 

Pl  =  R  g   T  (55) 

gel 

with  the  units  of  slugs  per  cubic  foot. 

From  the  survey,  the  velocities  were  known  as  a  function  of  radius, 

V(r).   Thus,  a  flow  rate  was  calculated  at  the  inlet  plane  by  numerical 

integration  of 


Wj  =  2  n  Oj  J     V(r)  r  dr  (56) 
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assuming  axisymmetr ic  flow  and  no  blockage  from  support  strut  wakes. 
This  flow  rate  was  found  to  be  three  to  four  per  cent  higher  than  the 
flow  rate  obtained  from  the  flow  nozzle  measurements.   This  comparison 
was  considered  satisfactory  for  the  range  of  measuring  capability. 

E.   REYNOLDS  NUMBER 

The  characteristic  length  for  the  Reynolds  number  (Re)  was  taken 
as  the  hydraulic  radius  of  the  inlet  annulus,  AR  ,  as  in  Ref.  2,  or 

2  (flow  area) 

^i  =  /  71   A : -r~T-  (57) 

1    (wetted  perimeter)  x   ' 


Therefore  , 


and,  substituting 


Dl  Vl  ^1 
Re  =  -J— ^ L-  (58) 


*  =  Pl  Vl  Al  =  Pl  Vl  n  (RT  "  4  }  (59) 


the  Reynolds  number  (Re)  can  be  expressed  in  dimensionally  correct  units 
as 

g  m  tt  (RT  +  Rjj  ) 

-7  2 

where  viscosity,  (j,  has  a  representative  value  of  3.9  x  10   lb  -sec/ft 

for  the  operating  range  of  85  to  95  degrees  Fahrenheit.   The  flow  rate, 

w,  is  in  lb  /sec,  and  the  radii  are  expressed  in  inches, 
m 
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V.   ANALYSIS  OF  TEST  DATA  FROM  REFERENCE  2 

Sovran  and  Klomp  [2]  made  a  study  of  more  than  one  hundred  annular 
diffuser  configurations,  and  their  published  data  was  of  a  form  suit- 
able to  be  analyzed  for  the  present  purpose.   Of  the  many  other  works 
published,  none  furnished  a  sufficiently  wide  cross-section  of  data 
which  could  be  readily  applied.   Data  from  Ref.  2  were  reduced  by  the 
Fortran  IV  program  shown  in  Table  Dl.   The  results  are  listed  in  Table 
D2  and  plotted  on  Fig.  16  for  a  Reynolds  number  of  approximately  6  x  10 
Each  set  of  data  points  represents  a  family  of  straight  wall  diffusers 
having  the  same  inlet  area  and  wall  divergence  angles  as  shown  in 
Fig.  16.   The  data  points  fall  within  a  band  on  the  plot  of  Q   cf  versus 
Q,  which  is  bounded  by  lines  representing  a  dimensionless  length 
description  of  the  diffuser,  L/AR.  An  investigation  of  the  values  of 
L/AR  for  existing  diffuser  designs  indicates  that  the  boundaries  shown 
on  Fig.  16  are  in  the  range  of  practical  designs.   The  points  falling 
below  L/AR  =  1.5  represent  designs  which  become  so  short  that  diffusion 
is  degraded,  and  those  above  L/AR  =  13  become  excessively  long  for 
practical  use.   It  can  be  seen,  that  as  Q  falls  below  approximately 
five,  there  is  an  abrupt  rise  in  the  value  of  fi  c,.,  which  is  repre- 
sentative of  the  entropy  rise  and  thus  of  an  increase  of  the  losses  in 
the  diffuser.   The  trend  of  the  curve  for  larger  values  of  Q  was  found 
by  considering  that,  as  a  diffuser  with  small  divergence  angles  becomes 
long,  its  characteristics  approach  t.;at  of  a  pipe.   Reference  9  gives 


the  value  of  the  coefficient  of  friction  for  an  annular  pipe  as 

f 


c     =  0.01  at  Re  =  6  x   10    .      A   pipe  with   L/AR  =   13  has  a  value   of 
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5  10 

£1    (SHAPE     PARAMETER) 


FIGURE    16 


Q  cf  vs.    Q  for  Reference  2   data, 
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Q  =  26  and  Q  c   =  0.26,  which  is  a  higher  value  of  Q  c   than  that 
occurring  from  diffuser  tests  at  Q   =  11.5  and  L/AR  =  13,  suggesting 
an  upward  trend  in  the  curve  for  high  values  of  fi.   This  trend 
indicates  that  a  minimum  of  Q  c   occurs  between  about  0=5  and 
Q  =  15,  possibly  at  a  value  of  Q  =  10,  for  an  L/AR  =  13. 

Vavra  [3]  calculated  values  of  0  c  ,  for  two-dimensional  subsonic 
diff users,  which  lie  within  the  same  band  as  the  data  plotted  on  Fig. 
16. 
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VI.   DISCUSSION  OF  TEST  DATA 

The  test  results  obtained  from  diffuser  Model  (1)  and  Model  (2)  are 
given  in  Table  II.   The  quantities  fi  c   for  Model  (1)  are  plotted  on 
Fig.  17  versus  Q.   Tests  were  made  at  different  flow  rates  to  investi- 
gate the  effect  of  Reynolds  number  on  the  location  of  points  in  the 
0  c,  vs.  Q  plot  for  the  different  configurations  of  Model  (1).   Figure 
17  shows  that  as  Reynolds  number  increases  the  quantity  Q  c   decreases, 
and  vice  versa.   This  effect  could  be  anticipated  from  the  trend  of  c 
against  Re  for  flows  in  annular  pipes,  as  shown  in  Ref.  9. 

The  curve  for  L/AR  =  6  and  the  curve  for  the  family  of  diffusers 

with  §   =  $ .  =  15  ,  from  Fig.  16,  and  the  curves  from  Fig.  17  were 

transposed  to  Fig.  18.   The  shape  of  the  curves  from  data  of  Model  (1) 

approximates  that  of  the  similar  family  of  diffusers  from  Fig.  16  for 

$  =  § .  =  15  .   However,  the  values  obtained  from  the  tests  of  Model  (1) 
o     1 

are  considerably  lower  than  those  from  Ref.  2.   This  discrepancy  may  be 

explained  by  the  fact  that  the  models  of  Ref.  2  were  constructed  of  wood 

and,  therefore,  had  a  higher  wall  roughness  than  the  cast  resin  models 

used  in  this  study.   The  value  of  c   for  a  wood  surface  would  be 

larger  than  that  for  a  resin  surface  with  lower  wall  roughness,  similar 

to  what  can  be  seen  from  the  Moody  diagram  for  pipes  of  differing 

roughness . 

The  values  of  Q  c   for  configuration  four  of  Model  (1)  became 

slightly  negative.   Values  of  Q  c   from  Ref.  2  data  were  also  negative 

when  C   became  greater  than  C   ..   For  an  adiabatic  process,  the  value 
pr  pri  v 

Q  c,  cannot  be  negative,  since  C   is  less  than  C   .  by  definition.   An 
f  pr  pri  J 

investigation  of  the  pressure  measurements  for  configuration  four  of 
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TABLE  II 
Model  (1)  and  Model  (2)  test  data 


MODEL 

CONFIGURATION 

fi 

Q  cf 

'P 

R  x  105 
e 

1 

1 

11.521 

0.0614 

.8793 

1.165 

0.0667 

.8837 

.879 

0.0566 

.8709 

1.486 

0.0632 

.8628 

1.514 

0.0573 

.8696 

1.518 

0.0569 

.8699 

1.510 

0.0610 

.8794 

1.201 

0.0560 

.8704 

1.526 

0.0859 

.8647 

.768 

2 

9.689 

0.0338 

.8677 

1.517 

0.0356 

.8683 

1.506 

0.0418 

.8734 

1.234 

0.0486 

.8808 

.775 

0.0529 

.8740 

.762 

3 

8.012 

0.0288 

.8699 

.766 

0.0135 

.8633 

1.509 

0.0170 

.8703 

1.208 

0.0200 

.8780 

.980 

0.0146 

.8587 

1.493 

0.0137 

.8607 

1.495 

4 

6.470 

-0.0162 

.8589 

1.502 

-0.0027 

.8756 

.755 

-0.0105 

.8772 

1.034 

-0.0128 

.8677 

1.265 

-0.0151 

.8589 

1.494 
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TABLE   II,    continued 


MODEL 

CONFIGURATION 

fi 

Ocf 

'p 

R     x    105 
e 

2 

1 

5.701 

0.0949 

.8226 

1.467 

0.1009 

.8279 

1.143 

0.1036 

.8362 

.741 

0.0941 

.8243 

1.466 
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Model  (1)  showed  a  fluctuation  of  the  instrument  readings  of  sufficient 

magnitude  for  C   to  be  larger  than  C 

pr  pri 

It  can  be  anticipated  that  the  curves  of  Fig.  16  will  be  shifted 
vertically  up  or  down  with  Reynolds  number  similar  to  the  condition  in 
Fig.  17.   The  values  of  0  for  the  point  of  sharply  increasing  Q  cf 
would  not  be  affected  by  a  change  in  Re,  since  the  entire  band  of  data 
points  would  be  displaced  in  the  vertical  direction  only.   To  avoid 
large  values  of  Q  c  ,  a  minimum  value  of  Q  =  5  should  be  used  for  most 
practical  designs. 

To  demonstrate  the  usefulness  of  the  design  parameter  Q,  for  com- 
plicated shapes,  Model  (2)  was  designed  as  outlined  in  Appendix  B  for 
a  value  of  Q  =  5.7,  and  L/AR  =  6.85.   Values  of  Q  =  6  and  L/AR  =  6 
were  chosen  for  the  original  design  point,  to  obtain  an  Q  cf  =  0.1. 

By  extrapolating  from  the  curve  for  L/AR  =  6  of  Fig.  16  to  the 
data  of  Fig.  17  to  show  the  effect  of  Re  and  the  difference  in  the 
model  surfaces,  the  heavy  dashed  curve  of  Fig.  18  was  obtained.   The 
data  of  Model  (2),  plotted  on  Fig.  18,  verified  the  location  of  the 
L/AR  curve  for  an  arbitrary  shaped  diffuser. 

The  efficiencies  for  both  Model  (1)  and  Model  (2),  plotted  on  Fig. 
19,  remained  essentially  constant  through  the  Re  range  of  this  study. 

Static  pressure  (p  )  on  the  inner  and  outer  walls  of  Model  (1)  and 
Model  (2)  was  non-dimensionalized  by 

C  -  ?tl  "  ^  (61) 

and  plotted  against  non-dimensionalized,  meridional  wall  length, 

X   =  1/L  ,  on  Figures  20  and  21,  where  1  is  the  meridional  distance 
w       °  ' 

along  the  wall  from  the  inlet  to  the  pressure  tap,  and  L  is  the 
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meridional  wall  length.   Values  for  C ,  X  ,  and  p  are  contained  in 
Table  III,  for  the  outer  wall  of  Model  (1)  and  for  both  the  inner  and 
outer  walls  of  Model  (2).   The  static  pressure  became  higher  than  the 
ambient  pressure  at  about  X    =  0.5  on  the  inner  wall  of  Model  (2), 
which  gives  a  Q    less  than  that  at  the  exit.   The  value  of  £  reaches  a 
minimum  between  X   =   0.7  and  X=0.9  and  increases  again  toward  C  =  0.293 
at  the  exit  plane.   This  condition  may  be  due  to  a  separation  of  the 
flow,  and  indicates  the  need  for  proper  wall  contouring  in  addition  to 
satisfying  the  criterion  of  Q. 
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TABLE  III 
Meridional,  static  pressure  values  for  Model  (1)  and  Model  (2) 


CONTOUR 

STATION 

X* 

p      (in.   water) 
gage 

c** 

OUTER 

1 

0.0 

- 13 .  54 

1.0 

CONFIG.    1 

Model    (1) 

2 

0.133 

-   8.79 

0.729 

3 

0.267 

-   5.46 

0.540 

4 

0.400 

-  3.50 

0.429 

5 

0.533 

-  2.24 

0.357 

6 

0.667 

-    1.45 

0.312 

7 

0.800 

-   0.82 

0.276 

8 

0.933 

-  0.34 

0.24^ 

9 

1.000 

0.0 

0.230 

OUTER 

1 

0.0 

-   9.66 

1.0 

CONFIG.    4 

Model    (1) 

2 

0.4 

-  4.90 

0.720 

3 

0.8 

-    1.45 

0.520 

4 

1.0 

0.0 

0.435 

*  X  -  The  dimensionless  length  1/L 


w 


**£  -  Dimensionless  pressure 
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TABLE  III,  continued 


CONTOUR 

STATION 

X 

p      (in.   water) 
w 

C 

INNER 

2 

0.107 

-7.97 

0.787 

Model    (2) 

3 

0.222 

-4.27 

0.558 

4 

0.336 

-2.14 

0.426 

5 

0.460 

-0.61 

0.331 

6 

0.578 

+0.21 

0.280 

7 

0.697 

+0.73 

0.248 

8 

0.817 

+  1.01 

0.231 

9 

0.933 

+0.85 

0.241 

10 

1.000 

0.00 

0.293 

OUTER 

1 

0.0 

-11.41 

1.000 

Model    (2) 

2 

0.121 

-10.91 

0.969 

3 

0.241 

-7.90 

0.783 

4 

0.366 

-5.18 

0.614 

5 

0.484 

-3.50 

0.510 

6 

0.605 

-2.40 

0.442 

7 

0.718 

-1.58 

0.391 

8 

0.823 

-1.15 

0.364 

9 

0.938 

-0.37 

0.316 

10 

1.000 

0.00 

0.293 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  paper  presents  a  design  parameter  Q,  for  diffuser  designs, 
which  takes  into  account  losses  due  to  the  geometry  of  the  diffuser. 
It  is  shown  from  reference  data  and  model  tests  that  there  is  an 
abrupt  increase  in  losses  for  diffusers  designed  for  a  value  of  Q 
less  than  about  five,  and  that  this  effect  is  independent  of  Reynolds 
number. 

The  parameter  Q  is  defined  for  incompressible,  adiabatic  flows.   It 
is  recommended  that  further  experiments  be  made  to  extend  the  data  for 
incompressible  flows  to  higher  values  of  fi,  and  that  an  investigation 
of  the  compressible  case  be  made. 

It  is  suggested  that  Model  (2)  be  improved  by  recontouring  the 
inner  wall  to  avoid  separation.   It  is  recommended  also  that  the 
present  inner  contour  be  cut  off  at  the  station  where  the  value  of  £ 
is  a  minimum.   Addition  of  static  pressure  taps  on  the  inner  wall  at 
the  inlet  plane  and  a  method  for  surveying  the  exit  plane  would  pro- 
vide data  to  better  describe  the  flow. 

It  is  recommended  that  the  method  of  measuring  pressure  be  investi- 
gated and  improved  to  give  more  precise  results  for  small,  fluctuating 
pressure  measurements,  prior  to  conducting  further  tests  on  diffusers 
with  Q  c  approaching  zero. 
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APPENDIX  A 

DETERMINATION  OF  DIFFUSER  EFFICIENCY 
The  diffuser  efficiency  was  determined  by  the  method  of  Vavra  [10] 
for  an  adiabatic  compression  process  with  friction.   The  polytropic 
efficiency  is  defined  for  the  differential  compression  from  p  to 
p  /  Ap  of  Figure  3  as 

dT. 

ip  s  "dF  (A1) 

The  dimensionless  flow  function  $  and  the  relation  for  the  poly- 

c  r       j 

tropic  exponent,  n,  from  Ref.  10  are 

-y  ,     I.    .   2/n      ^ 


,-^-L=^[^(^'"-(|)n]        m 


$   =  w  -&  T. 
c     \gQ      1   PjA2   J  7- 


and 

7  r\ 


2_ 


1-7  d-Tip) 


(A3) 


where  w,  T   ,  T1 ,  p_  and  p,  are  either  measured  or  computed  for  the 
particular  diffuser  design. 

For  an  incremental  change  in  pressure  from  p.  to  p_ ,  the  pressure 
Py  may  be  expressed  as  p1  /  Ap,  where  Ap  is  small  with  respect  to  p, . 
The  pressure  ratios  in  (A2)  can  be  expressed  as 


*-0  ♦£) 


El 

pi 


By  the  binominal  expansion  theorem, 

2/n  2/n 

(*T)    _(l  +  iE)    -l  +  »*B+...  (A5) 

M^      V    pj^         n  Pl 
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and   similarly 


n+1 
Pox  n 


© 


n      p 


(A6) 


and  by  assuming  Ap  to  be  small,  the  (Ap/p  )   and  higher  order  terms  are 
much  smaller  than  Ap/p  and  are  omitted. 
From  (A2) 


n+1 


T  2/n 

»  2  2=1  =  Ji  (?£)      _  Cp-T)   n 

c   27      tx    \P]J        vPly 


(A7) 


Substituting  expressions  (A5)  and  (A6) ,  and  solving  for  n,  gives 

T 


^2 
Pl 


2  =—  -   1 

-     Tl 


n!=ZziS2.Tti+i+An 


ly       c 


(A8) 


Equating  (A3)  and  (A8)  gives 

T 

Ap. 


T|   = 
P 


'1 


L2T 


?-] 


&(£-0(?£+0- 


(A9) 


with 


i  -  1  +  ^  M? 


2    1 


(A10) 


the  polytropic  efficiency  is 


Ap_ 
pl  L 


1  + 


(7-1)  M^] 


$   -  >M,   2  — c  +  1 1 
c   'Up.    J 


(AH) 


62 


APPENDIX  B 

DESIGN  OF  DIFFUSER  MODEL  (2) 

From  the  results  of  Ref.  2  and  the  test  data  of  Model  (1),  a  design 
point  of  Q  =  6  and  L/AR  =  6  was  chosen  from  Fig.  18  for  Model  (2).  The 
inlet  area  of  Model  (2)  was  the  same  as  for  Model  (1). 

From  equation  (42)  of  the  main  text, 

C   .  -  C   =  Q   c 
pri    pr      f 


or,    from  equation    (40),   and   a   value   of  fi  c      =  0.1, 


pr 


(Bl) 


For  the  range  of  C   from  0.65  to  0.72  for  a  L/AR  =  6,  equation  (Bl) 

Pr 

gives 

A2=-l-T70+C   )  <B2) 

pr 

2 

giving  a  range  of  A_  from  39  to  51  in  . 

An  approximate  radius  was  chosen  to  give  a  reasonable  axial  opening 

2 

at  a  radial  exit  plane  for  A„  =  45  in  . 

Since  the  outer  surface,  E,  of  Fig.  22  was  most  critical  from  the 
standpoint  of  flow  separation,  it  was  designed  first.   A  smooth  curve 
of  the  desired  shape  was  sketched  in  for  surface  E.   To  approximate 
the  curve  sketched,  a  series  of  circular  arcs  were  constructed  which 
provide  for  simple  tooling  layouts. 

Figure  23  shows  the  method  by  which  two  circular  arcs  are  fitted  to 
satisfy  the  condition  that  both  have  centers  on  the  same  line  at  their 
tangent  point  (4),  and  that  they  fit  into  a  given  length,  L. 
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FIGURE  22 
Design  of  Model  (". ) 
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FIGURE  23 
Fitting  of  circular  arcs 
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The  largest  radius,  R.  ,  is  chosen  to  give  the  desired  effect  of 
the  original  sketch.   A  curve  is  then  drawn  from  point  (1)  to  point 
(2),  where  (2)  is  an  approximate  point  of  the  desired  curvature  change 
Point  (2)  defines  a  and  b  for  any  given  distance  d.   The  center  of  R- 
becomes  point  (3). 

From  geometry,  R~  and  R'   can  be  defined  as 


R9  =  j  9n  (B3) 

|l  -  <a/R.r 


and 


ab 
R'2  =  L  -  a  /  r— - j  (B4) 

K  "  a 


Point  (2)  may  be  iterated  until  R.  =  R'9-   The  centers  of  R1  and  R„ 
will  then  be  on  the  same  radius  at  their  tangent  point  and  the  curve 
will  fall  within  the  axial  length  L  and  height  d. 

The  same  procedure  as  described  above  was  extended  to  surfaces  of 
more  than  two  circular  arc  segments,  such  as  surface  E,  of  Fig.  19,  by 
locating  the  arc  centers  E-l,  E-2,  and  E-3 .   Surface  F  was  constructed 
in  a  similar  manner. 

With  areas  A  and  A„  known,  a  linear  distribution  of  flow  area  from 
station  (1)  to  station  (2)  was  chosen.   Circles  of  the  diameter  which 
gave  the  desired  area  at  a  given  meridional  station  were  then  con- 
structed along  an  approximate  mean  streamline  as  shown  by  the  example 
circles  on  Fig.  22.   Surface  B  was  then  sketched  in,  tangent  to  the 
area  circles,  and  fitted  with  circular  arcs  as  described  above.   The 
flow  channel  was  then  checked  for  uniformly  increasing  area  as  a  function 
of  length,  and  adjusted  to  prevent  any  convergence  in  the  flow  area. 
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APPENDIX  C 

COMPUTER  PROGRAM  OMEGA 
Program  "Omega"  is  a  modification  of  program  "Diffuser"  used  by 
Beck  [6].   The  program  is  written  in  Fortran  IV  for  the  IBM  360-67 
digital  computer.   The  program  was  established  to  compute  the  values 
of  flow  rate,  Re,  Q,  Q  c  and  r\     by  the  methods  described  in  the  main 
text.   Table  CI  is  a  listing  of  the  program. 

Cl.   Assembly  of  Input  Data 

The  field  specifications  for  the  input  data  are  shown  in  Fig.  24. 
Card  one  of  the  data  deck  is  a  72  space  "Hollerith"  statement  for  data 
identification.   The  content  of  this  card  is  printed  immediately  below 
the  output  title.   Card  two  is  a  two  digit  I-format  number  indicating 
the  number  of  data  points  in  the  deck.   Card  three  is  the  input  of 
inlet  and  exit  diameters  and  the  value  of  fi.   Cards  four  and  five  con- 
tain the  data  recorded  for  one  data  point  and  must  be  repeated  for 
each  point.   No  other  cards  need  be  repeated. 

All  floating  point  numbers  may  be  placed  anywhere  in  their  field. 
I-format  numbers  must  be  right  justified. 

The  definitions  of  the  program  input  variables  are  given  below. 

Fortran 

Symbol         Definition  Format 

NUM  Number  of  sets-  of  data  points  I  2 

D01  Diameter  of  outer  wall  station  one  F  10.0 

DI1  Diameter  of  inner  wall  station  one  F  10.0 

D02  Diameter  of  outer  wall  station  two  F  10.0 
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Fortran 

Symbol 

D02 
DI2 

OMEGA 

IRUN 

MO 

IDAY 

PATM 

PINOZ 

DPNOZ 

TNOZ 

TEMP 

PI 

Ql 
TT1 

TRM 


Definition 

Diameter  of  outer  wall,  station  two 

Diameter  of  inner  wall,  station  two 

Value  of  shape  parameter 

Number  of  run 

Month 

Day 

Atmospheric  pressure 

Static  pressure  upstream  of  nozzle 

Pressure  differential  across  nozzle 

Total  temperature  upstream  of  nozzle 

Temperature  read  at  barometer 

Static  pressure  at  station  one 

Dynamic  head  at  station  one 

Total  temperature  measured  upstream 
of  station  one 

Ambient  room  temperature  at  exit 


Format 

F 

10.0 

F 

10.0 

F 

10.0 

I 

3 

I 

3 

I 

3 

F 

10.0 

F 

10.0 

F 

10.0 

F 

10.0 

F 

10.0 

F 

10.0 

F 

10.0 

F 

10.0 

F  10.0 


C2.   Output 

The  program  output  is  in  three  sections.   The  first  section  is  a 
printout  of  data  as  read  from  the  input  cards  as  a  check.   The  second 
section  is  the  calculation  of  flow  rate  which  is  divided  into  input 
values  converted  by  the  program  to  dimensionally  correct  units,  and 
output  of  flow  rates  computed  by  the  program.   The  third  section  is 
also  divided  into  input  values  converted  by  the  program  and  output 
values  of  interest  in  evaluation  of  diffuser  performance. 
Table  C2  is  a  sample  output  of  the  program. 
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CARD 

1 

72 

72   Space  Hollerith   Statement 

Data   Identification 

Blank 

CARD  2 


1  2 


I  3 


NUM 


Blank 


CARD  3 


1               10  !  11              20 

21              30 

31              40 

41              50 

F10.0 

F10.0 

F10.0 

F10.0 

F10.0 

D01 

DI1 

D02 

DI2 

OMEGA 

Blank 

CARD  4 


1       3 

4         6 

7        9 

10        19 

20       29 

30       39 

40       49 

50        59 

I  3 

I  3 

I  3 

F10.0 

F10.0 

F10.0 

F10.0 

F10.0 

IRUN 

MO 

IDAY 

PATM 

PINOZ 

DPNOZ 

TNOZ 

TEMP 

Blank 

CARD  5 


1            10 

11            20 

21           30 

31            40 

F10.0 

F10.0 

F10.0 

F10.0 

PI 

01 

TT1 

TRM 

Blank 

FIGURE  24 
Input  field  specifications  for  Program  Omega 
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APPENDIX  D 


REDUCTION  OF  DATA  FROM  SOVRAN  AND  KLQMP,  REFERENCE  2 
Data  of  Ref.  2  were  reduced  by  a  Fortran  IV  program  to  a  form 
useable  for  this  study.   Table  Dl  is  a  listing  of  the  program.   The 
diffuser  description  of  Ref.  2  was  the  same  as  in  Fig.  10. 
The  definition  of  the  input  variables  are 


Fortran  Symbol 

ALFAH 

ALFAT 

DIMLL 

ARATIO 

RRATIO 

RT1 

CPR 


Description 

The  divergence  angle,  cl,,  of  the  inner  wall 

The  divergence  angle,  (T,  of  the  outer  wall 

The  dimensionless  length,  L/AR. 

Area  ratio,  A  /A 


Radius  ratio,  R^/R^ 
Radius  of  outer  wall,  station  one 
Coefficient  of  pressure  recovery 
The  axial  length  was  computed  from  the  geometry  and  a  numerical 

integration  was  made  to  obtain  Q.   The  ideal  pressure  recovery 

coefficient  was  computed  from  the  area  ratio  and  used  to  compute  Q  c 

The  results  are  contained  in  Table  D2 ,  where 

Fortran  Symbol  Description 

OMEGA  The  shape  parameter,  fi 

CFBAR  Average  local  coefficient  of  friction,  c 

DIMLL  The  dimensionless  length,  L/AR 

OCF  The  product,  Qc 
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TABLE  D2 
SOVRAN    &    KLOMP,    GMR-511,     DATA 
OMEGA  CFBAR  DTMLL  OCT 


3o  64  54 

0o0027 

2,2800 

0o0097 

5o2663 

0o0054 

3,6300 

0,0283 

7o2076 

0o0117 

5o6600 

0o0341 

9Q1634 

OoOl 25 

8,3700 

0oll46 

9o9462 

Oo0136 

9,7?00 

0ol351 

10o63?l 

0o0137 

llo0700 

0,1461 

llo5207 

Co0132 

3.3,1000 

0ol519 

3G6454 

Co0013 

2o28C0 

0o0047 

7o2076 

0o0119 

5„6600 

O0O86I 

702C"76 

Co  CI  24 

5,6600 

Co0891 

7o2C76 

0o0143 

5o660C 

Col031 

9ol634 

Co  01 35 

8,3700 

0ol236 

9ol634 

0oC143 

8o2700 

0ol30fc 

1U5207 

0o0147 

13,1000 

Col68c 

2o96?9 

0o0221 

2o2  600 

0,0655 

3o4199 

0o0245 

2o93C0 

O0II8I 

3o7999 

0o0291 

3o6100 

0oll05 

4048  38 

Co036° 

4o96CC 

0,1656 

4o5598 

0o0291 

5o6t00 

Ool32  8 

4o63  5  8 

0o0406 

6o3200 

0ol88^ 

5o0918 

0o0391 

8o3400 

Col993 

5o3198 

0o0372 

10,3700 

Col978 

505478 

0oC388 

13,C80C 

0o2151 

2o9639 

CoC262 

2o260C 

0o0775 

3o7999 

0o0306 

3o6l00 

Go  11 65 

3o7999 

0o0391 

3o6100 

0ol435 

3o7999 

0o0375 

3o61C0 

0ol425 

3o7999 

Oo0288 

3o6100 

0ol475 

3o7999 

0o0430 

3o6100 

0ol635 

405598 

0o0296 

5,6400 

0ol348 

4o5598 

0o0344 

5o6400 

0ol568 

4o5598 

0o0368 

5o6400 

0ol678 

4o5598 

Co0421 

5o640C 

Co  191 8 

2o71^2 

0o0511 

2o25C0 

C1387 

3ol363 

0o0575 

2o9300 

Co  1804 

3o4570 

0o0631 

3o6000 

Co2182 

2o42^0 

0oC957 

2o2400 

0o?320 

2o7^2C 

Ool057 

2o9100 

0,2900 

2o9811 

0ol098 

3o5  900 

Co  3 2  74 

3o4800 

Co  1070 

6„2  90  0 

0o3722 

3o7C44 

0ol0^6 

9o6800 

0o3875 

3o7914 

Co0983 

13o0600 

Co3728 

2o7884 

0o0518 

2o3600 

0ol444 

2o9829 

0oC547 

2o7100 

Co  1633 

3o0478 

0o0570 

3o05CC 

Co  173  6 

3o3071 

0o05Sl 

3o4000 

0ol856 

3o5017 

0o0556 

3o7500 

Co  1949 

3o6314 

0o0  514 

4,4400 

O0I866 

3o8908 

0o0471 

5ol400 

0,1833 

4ol501 

0o0427 

5o8300 

0,1772 

77 


TABLE  D2,    continued 

SOVRAN    £    KLOMP,    GMR-?11,     DATA 
OMEGA  CFBAR  DIMLL  OCF 


2o7884 

Oo058? 

2o^600 

0ol62^ 

2o9829 

Oo0624 

2o7100 

0ol863 

2o2947 

0ol408 

2o3  400 

0o3?31 

204358 

Ool429 

2,6800 

0o3480 

2o5578 

0o1468 

3o0300 

Go3756 

2o6602 

0ol48^ 

3o3800 

0,3949 

2o7467 

C  o  1  5 1  5 

3o7300 

Go  41 6 3 

2o88?0 

0o14  8'- 

4o420C 

0,4276 

209835 

0ol^4?. 

5ol2G0 

0o43C5 

3o0580 

0ol297 

5o80CC 

0o4273 

?ol696 

Col6°l 

2o3300 

Co3669 

2o2974 

0ol723 

2o6  800 

0,3959 

2o4C30 

0ol738 

3o0300 

0o4177 

2o^8  86 

0ol761 

3o3700 

0o4382 

2o5627 

0ol768 

3o7200 

0,45?0 

lo97^8 

Go  21 88 

2o32C0 

Co4321 

2o0752 

0o2247 

2o6600 

Co4662 

2ol5c3 

0o2278 

3o01C0 

0o4919 

2o2281 

0o2309 

3o3600 

0o5146 

2o2851 

0o2322 

3o7100 

0,5305 

2G4784 

0  o  2  3 1 0 

5o79G0 

0o5725 

2o8  5Q0 

0o0378 

2o42G0 

G0IO8I 

3o8475 

0o0298 

3,9200 

Go  1147 

4o7895 

Co0249 

60IIOO 

0oll90 

5o5813 

Go  02 10 

9o1700 

0oll72 

lo8581 

-0oCC3C 

1,4300 

-0,0C57 

2o4553 

Oo0153 

2o3C00 

0o0376 

3o0217 

0oC276 

3,6000 

0o0833 

3o4734 

0o0328 

5o3400 

Go  1139 

lo8581 

-0oGC36 

lo4300 

-C0C67 

lo8581 

0o0126 

lo430G 

0o0233 

204553 

0o0124 

2,3000 

0o0306 

3o0217 

0o0226 

3o60G0 

0o0683 

3o4734 

0oG230 

5o34GC 

00079^ 

304734 

0o0269 

5,3  400 

Go  1280 

2o3698 

-GoC155 

1,^700 

-Go0368 

3o3687 

-0oCG57 

2o3400 

-0o0191 

4o4Q30 

-Go  GO  12 

3,6^00 

-0o0053 

5o53A8 

0o0037 

5,3800 

0o02G6 

2o3698 

-Co0109 

lo4700 

-Oo0258 

5o5248 

0o0102 

5o3800 

0o0566 

lo6569 

-0o0004 

lo5000 

-0,0006 

lo8824 

OoOOOl 

lo9500 

0oG0C2 

2oC503 

0o0003 

2o4000 

0,0006 

2ol776 

OoCGOO 

2o8500 

OoOOGO 

2o2757 

-0o0002 

3,3  000 

-0o00C4 

2o3525 

OoOGOl 

3o7500 

OoOOC? 

2o4134 

0ol735 

4o2  000 

0,4186 

2o4631 

0ol723 

4o6  6G0 

Go  4244 

2o5C24 

0ol725 

5oll00 

Co4316 
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TABLE  D2,    continued 

SOVRAN    £    KLOMP,    GMR-511,    DATA 
OMEGA  CFBAR  DIMLL  OCF 


2o5346 

Ool732 

5o5600 

Co439C 

U6569 

0o1409 

lo5000 

0*233*- 

2o5346 

0o1846 

5o5  600 

0o4680 

2o8?.54 

0oC021 

2o0100 

CO  0  5  8 

3o5666 

0o0045 

2o8600 

0o0162 

4o2207 

O0OIO8 

3o8100 

Go  045 7 

4o7081 

Oo0122 

4o7100 

0o057* 

5o4215 

Oo0135 

60  5  20  0 

0,>073  4 

5o6828 

Oo0135 

7o^200 

0o0768 

5o8977 

OoOl 20 

8o3100 

0o0705 

2o8154 

Oo0038 

2o0100 

0o0!C8 

2o8C53 

OoOC67 

2o0000 

O0OI88 

lo4279 

0o2348 

U560C 

0o3353 

lo5757 

0o3059 

2o0200 

Co4820 

lo6803 

0o3143 

2o4900 

0o5281 

lo7532 

Oo3205 

2o9500 

0o5  619 

lo80  75 

Oo3232 

3o4200 

0o5841 

lo84"72 

Oo3242 

3o8800 

0o5989 

U4279 

C«2313 

lo5600 

0o3303 

lo8^72 

Oo3291 

3,8800 

Co6079 

2o2047 

Ooll29 

2oOPOO 

Co2488 

2o4275 

0o0844 

2o5400 

0o2048 

2o85"79 

Oo0903 

3o9300 

0o2582 

2o2C47 

Ool201 

2oC80  0 

0o2648 

2o8579 

Oolll3 

3o9300 

0o3182 
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